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factor for IRAS is near unity. However, one should be cautious in claiming that these results provide 
definitive proof of = 1. Similar estimates of (3 using optical surveys have yielded lower values 
(0.25 < (i < 0.65) for the optical dipole [19] and (J3 = 0.37 ± 0.2, (3 = 0.5 ± 0.06) in comparisons of 
optical samples with peculiar velocities (E. Shaya and M. Hudson respectively, these proceedings)) 
suggesting that at least the ratio of the IRAS and optical bias factors may differ from unity. 

Although the estimate of (3 depends quite sensitively on the normalization of the power spectrum, 
the method gives a surprisingly robust measurement of the shape of the power spectrum with T ~ 
0.16 ±0.05. This result is somewhat lower that the value quoted by Feldman, Kaiser, & Peacock [7] of 
r = 0.31 ± 0.08 in an analysis of the QDOT power spectrum and that of Mo, Peacock and Xia [18] of 
T = 0.32 ± 0.07 deduced from the cluster-galaxy cross-correlation function. These measurements all 
argue strongly for a real space power spectrum that is significantly steeper on scales £ 30 h~ 1 Mpc 
than the canonical CDM model of T = 0.5 (ruled out in all three analyses at £ 3-<7 level!). In the 
usual parlance, this can be phrased as the standard CDM model having an insufficient ratio of large 
to small scale power. Models containing a combination of hot and cold dark matter may provide one 
acceptable revision of the CDM model (cf. Pogoysan, these proceedings and references therein) as 
may the introduction of a cosmological constant [6]. 

The SHA analysis of redshift distortion presented here is well suited to near full-sky catalogs. The 
radial dependence of the selection function in flux-limited catalogs leads to a very compact expression 
for the distortion in terms of the harmonic power spectrum. The technique outlined in this paper, 
however, rests on the choice of an unspecfied radial weighting function, which induces the distortion. 
At first glance this arbitrariness may seem unappealing, however, it does allow the flexibility to 
optimize scheme for maximum distortion on a given scale. Alternatively, one can follow the approach 
used here and simply use a set of multiple weighting schemes which probe, albeit not optimally, a 
variety of different scales. Perhaps a more natural way to formulate the distortion is in terms of an 
expansion of the redshift space density field in spherical harmonics and orthogonal radial functions 
(see 0. Lahav contribution to these proceedings). Spherical coordinates have the advantage that 
the redshift distortion couples only the radial modes thereby offering the possibility of removing the 
distortion using regularized inversion techniques. This work is currently in preparation and will be 
presented at a later date. 

Acknowledgements. I would like to thank my collaborators Ofer Lahav and Caleb Scharf as well 
as Michael Strauss, John Huchra, Marc Davis, and Amos Yahil for permission to use the 1.2 Jy IRAS 
survey data prior to its publication. 
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Figure 1: Likelihood contours derived from the SHA of the 1.2 Jy IRAS redshift survey. In this 
plot, the normalization of the power spectrum (<jg) has been held fixed at 0.69 (the value determined 
independently from the real space correlation function). The distortion is therefore a sole function of 
(3 and the shape of the power spectrum (parametrized by T). The contour levels correspond to the 
appropriate A% 2 values for the confidence levels labeled in the figure for a likelihood function with 
two degrees of freedom. 
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the same approach used by Yahil et al. [27] in their reconstruction of the IRAS peculiar velocity field. 
Moreover, regularized mask inversion of the harmonics shows the effect of incomplete sky coverage 
is small for the IRAS 1.2 Jy sample geometry and / < 10 (Lahav et al. 1993, in preparation). The 
selection function, (f>(r), is computed using the techniques of Yahil et al. [27]. In practice the selection 
function is computed by assigning luminosities based on galaxy redshifts, not distances; fortunately, the 
shape of the (f>(r) is insensitive to the flow model used to correct for the effects of peculiar velocities 
[21]. Moreover, Monte-Carlo simulations have shown the results to be unaffected by the expected 
measurement redshift errors in the survey [11]. 

Since the distortion term in the harmonics given in Equation (9) is proportional to both the 
amplitude of the power spectrum (<tJ) and (3, we expect a strong covariance in the likelihood estimates 
of these two parameters. Fortunately, the normalization of the real space power spectrum is fairly well 
established from analyses of the real space correlation function [10] and projected spherical harmonics 
[22]. As a first attempt, we computed the maximum likelihood estimates of (3 and T (the shape 
parameter of P R (k)) holding the normalization of P R (k) fixed at the value determined by Fisher et 
al. [10], us = 0.69 ± 0.04. The derived maximum likelihood values in the case are (3 = 0.94 ± 0.17 and 
r = 0.17 ± 0.05 and the corresponding likelihood contours are shown in Figure 1. 

Allowing (Ts to be a free parameter increases the statistical uncertainty in (3 and leads to the 
following maximum likelihood values: (3 = 0.47 ± 0.25, T = 0.15 ± 0.05, and a 8 = 0.81 ± 0.06. Table 2 
lists the correlation coefficients between the parameters. Table 2 shows the very strong anti- correlation 
between ag and (3. The shape parameter T, unlike ag is positively correlated with (3. Since the overall 
amplitude of the distortion is fixed by the data, an increase the predicted amount of large scale power 
(decreasing T) results in a lower inferred value of (3. The recovered value of ag when it is treated as a 
free parameter is just within the 3 — a error limits of the ag determined from the correlation function 
although its statistical error ag is 50% larger. Given the strong degree of (3-ag covariance, the more 
reliable value of (3 is probably the one derived by treating ag as a fixed parameter set by the real space 
correlation function (ag = 0.69). 

6 Discussion 

The value of (3 derived from the 1.2 Jy sample using the SHA technique is dominated by uncertainties 
in the amplitude of the power spectrum. The value of (3 = 0.94 ± 0.17 determined by holding ag fixed 
at 0.69, is in agreement with previous estimates from IRAS redshift surveys. Dekel et al. [4] found 
(3 = 1.28^q'5q (2ct) by comparing the density field of the 1.936 Jy IRAS survey [25] with the density 
field inferred from direct measurements of the peculiar velocity field using the POTENT algorithm 
([1], [3]). Kaiser et al. [16] found [3 = 0.9 jlo'f ° (la) in a comparison of the predicted IRAS peculiar 
velocity field with direct measurements using the deeper, albeit sparser, QDOT redshift survey. J. 
Roth (these proceedings) finds (3 ~ 0.85 ± 0.30 based a comparison of directly measured peculiar 
velocities with the reconstructed maps of the IRAS 1.2 Jy peculiar velocity field (derived using the 
analysis of Yahil et al. [27]). Measurements based on the predicted acceleration of the Local Group 
using IRAS catalogs have also yielded high estimates 0.4 < (3 < 1.0 ([20], [26]). Hamilton's [14] 
analysis of redshift distortions in the correlation function for the 1.936 Jy IRAS survey gave a slightly 
lower value, f3 = 0.66^:22 ( 

Ict), in good agreement with the value we obtained ((3 = 0.47 ± 0.25) by 
allowing ag to be a free parameter. 

These results provide evidence of a high, near closure, value of fi on large scales if the bias 
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covariance matrix, Aij, are £ 10% of the largest diagonal elements; thus the four spatial windows 
give roughly independent estimates of the redshift distortion. We will use these windows in all the 
analysis (both A-body and actual data) discussed below. 

4 Monte-Carlo Tests 

In order to check the validity of our formalism, we applied the likelihood formalism to the real and 
redshift space weighted harmonics computed from an A-body simulation of a standard Cold Dark 
Matter universe characterized by Cl h = 0.5. The simulations evolved 64 3 particles in a box of 180 
/i _1 Mpc using the P 3 M algorithm until rms variance of the density field in a sphere of 8 /i _1 Mpc 
reached <jg = 0.61; further details of the simulations can be found in Frenk et al. [12]. In an effort 
to mimic current observational data, we extracted nine mock galaxy catalogues designed to closely 
match the properties of the 1.2 Jy IRAS redshift survey (cf., [25], [8]). Galaxy candidates were chosen 
as unbiased tracers of the underlying particle distribution and therefore (3 = 1 in the mock catalogues. 
The procedure for extracting these catalogues is described in Gofski et al. [13] and Fisher et al. [9]. 

The A-body simulations provide a convenient way to test the formalism outline in the previous 
sections since both the amplitude and shape of the power spectrum are known a priori. Thus a simple 
test is to compute the weighted redshift harmonics for each of the mock IRAS catalogs and perform 
the likelihood analysis for the best value of (3 holding the shape and amplitude of the power spectrum 
fixed at their correct values. Since our analysis is valid only in the linear regime, we restricted the 
likelihood computation to / < 10. Table 1 shows the resulting estimates of (3 for the nine mock 
catalogs. The mean value of (3 for the nine mock samples is 0.94 ± 0.2 which is in excellent agreement 
with the actual value (3=1. Thus, when the correct power spectrum is used the likelihood estimate 
for (3 is unbiased; moreover the formal uncertainty in a single realization (given by AC = —1/2) is 
comparable to the scatter over the nine realizations, i.e., A/3 = 0.2. 

5 Application to the 1.2 Jy IRAS Survey 

We now proceed to apply the formalism described in the previous section to redshift survey of 5313 
IRAS galaxies flux- limited to 1.2 Jy at 60/um, selected from the IRAS database ([24], [8]). The 
sample covers 87.6% of the sky and is complete for \b\ > 5° with the exception of a small area of the 
sky not surveyed by IRAS. The SHA analysis discussed above relies on complete 4ir steradian sky 
coverage. Although statistical corrections can be applied to the harmonics with partial sky coverage 
([22], [23]), we have adopted a simpler method of dealing with incomplete sky coverage. We have 
interpolated the redshift data through the plane in a way which smoothly continues structure; this is 



We stress, however, that in linear theory the choice of reference frame (i.e., the value of V 6 s ) affects 
only the value of the dipole harmonic and that for / > 2 the results are independent of the assumed 
motion of the observer. The rms averaging procedure in Equation (9) treats the observer as a typical 
point in space and therefore assumes that the observer's peculiar motion is not unusually high or low. 
This point should be kept in mind if the motion of the Local Group is atypical. In practice, these 
complications can be avoided by working in the comoving frame of the observer, i.e., the reference 
frame in which \ bs vanishes or omitting the dipole harmonic in the analysis. 

The predicted rms redshift harmonics depend on three quantities: the value of (3 and the shape 
and normalization of the power spectrum. The shape of the power spectrum can be conveniently 
parametrized by a series of phenomological CDM models with with varying T = £l h (e.g., [5]). For 
the normalization, we adopt the standard <jg convention corresponding to the variance in spheres of 8 
/i _1 Mpc. 

3 Likelihood Estimator 

Suppose one computes a set of harmonics, a] m , for a set of different weighting functions, /W(r), for 
{i = 1, . . .N}. If the underlying density field is Gaussian, then the real and imaginary parts of the 
spherical harmonic coefficients are independent Gaussian variables. In this case, the likelihood (defined 
as the negative logarithm of the joint probability) is given by 

ln£[/3,cr 8 ,r] ee -^^^inldetA^I-^^M^o^W] 

I 

~ Y, A ^T,[Ma}JRe(aiJ + Im(a}JIm(a^ , (11) 

i,j m = l 

where Aij is the covariance matrix 

Aij = ( a im a Z)TH + {a) m a\* m )sN , (12) 

with 

oo 

<«LOth = lj dkk*p R {k) (vfhk) + ^fhk)) x 



(9*\k) + fMfV\kj)* (13) 

and 

oo 

(*} m <JsN = J drrV(r) [/W(r)/^(r)] , (14) 
o 

In Equation (11) Re(ai m ) and Im(ai m ) refer to the real and imaginary parts of a\ m [23]. The super- 
scripts on ^ R (k) and ^ c (k) refer to the window function of the corresponding weighting function. 
Since the harmonics at different / are independent, the total likelihood for / < l max is just given by the 
sum of the likelihoods in Equation (11) for each / up to l max . Even if the primordial density is a true 
Gaussian field, nonlinear evolution will lead to departures from Gaussianity on small scales. Conse- 
quently, we expect the formulation of the likelihood given above to be correct (under the assumption 
of Gaussian primordial fluctuations) only for low order harmonics which probe fluctuations still in the 
linear regime. 

The number and functional form of the weighting functions, f^\r), used to construct the har- 
monics is arbitrary. Clearly, one would like to pick those functions which simultaneously maximize 
the distortion while minimizing the noise. After experimenting with several different choices of weight 
functions we decided to use four Gaussian windows centered at 38, 58, 78, and 98 h~ 1 Mpc each 
with a dispersion of 8 h~ 1 Mpc. For this choice of weight functions, the off-diagonal elements of the 



There are several comments to be made regarding the validity of Equation (4). First, we have 
assumed that the summation in Equation (2) is carried out over all galaxies in a flux limited redshift 
survey. In this case, the integrals in Equation (4) extend over all space and there are no "surface" 
terms arising from the deformation of the boundary of the integration region that occurs in the 
transformation from redshift to real space. Second, there is an apparent absence of terms involving 
the derivatives of the selection function; these terms would be manifest if we, like Kaiser [15], took 
f(s) to be the special case of l/cj)(s). 

The expansion in Equation (4) can be simplified by expanding the radial peculiar velocity field 
dependence, U(r), in spherical harmonics: 

U ( r ) = £~2 £(*')* / d " k f jl(kr)Y lm (k)Y^(r) . (5) 

In Equation (5) and below, j'i(kr) = dji(kr) / 'd(kr) refers to the first derivative of the Bessel function. 
Equation (5) can be derived in a straightforward way by using the linearized continuity equation to 
relate the velocity to the density field, = — i(3S^k/k 2 , and the Rayleigh expansion of a plane wave 
in spherical waves. Using Equation (5), Equation (4) can be written (for / > 1) as 

aL = S7 d3kS ^?( k ) + ^f(MY lm (V , (6) 

where 

CO 

*f(A;) = J drr^f^kr) (7) 
o 

describes the real space contribution to the harmonics and 

oo 

¥f (*) = \ j drr 2 #r)^ (j>( kr) _ ^n) , (8) 

o 

is a "correction" term which embodies the redshift distortions. In Equation (8), 8n is a Kronecker 
delta which contributes only to the dipole (/ = 1) harmonic. In the case of equal weight, df/dr = 0, 
^f(k) vanishes and Equation (6) leads to the real space expression for the harmonics derived by Scharf 
et al. [22]. The 1/k scaling of the distortion window, ^f(k) shows that it will generally be much 
more sensitive to long wavelength perturbations than the real window ^f(k). 

The expected linear theory rms value of the harmonics is given by ensemble average of the square of 
Equation (6). Using the statistical independence of Fourier waves 2 , {Sffi^T) = (2'K) 3 P R (k)S^ 3 \'k — k'), 
this yields 

oo 

<l«fj 2 ) = - [dkk 2 P R (k)\^(k) + ^f(k)\ 2 . (9) 
o 

In real data the square of harmonics in Equation (9) will have a discreteness or "shot" noise contri- 
bution; this can be modeled by adding 

oo 

(\ai m \ 2 ) SN = j drr 2 cp(r)[f(r)] 2 (10) 
o 

to Equation (9). 

The / = 1 or dipole distortion has two contributions: one from the external dipole moment of 
the velocity field around the observer and one induced by the motion of observer itself. The latter 
distortion is caused by V 6 s ^ and is the origin of the so-called "rocket effect" (c.f., [15], [17], [26]). 



2 This assumes that the survey contains many independent modes of the wavenumber in question and usually phrased 
somewhat loosely as the "fair" sample hypothesis. 



respectively; we also will use R and S super and subscripts when referring to quantities in real and 
redshift space respectively. 

The quantity (3 naturally appears in Equation (1) since it is the controls the amplitude of the 
peculiar velocity field; in linear theory, (3 is the proportionality factor relating the velocity and density 
fields with V • v = —H (3S. The dependence of the clustering distortion on the cosmic density 
parameter has lead to a variety of methods (e.g., [14], [2] and references therein) aimed at exploiting 
the measurements of the redshift distortion as a way to determine fi . Here we present a technique of 
investigating redshift distortions based on a linear expansion of the galaxian density field in spherical 
harmonics. The radial nature of the peculiar velocity field leads to a simple and compact formulation 
of the redshift distortion in spherical coordinates. The work presented here is an extension of the 
method developed in Fisher, Scharf, & Lahav [11]. Spherical harmonic analyses of the density field 
(hereafter SHA) have been shown to be a successful way to quantify large scale structure and have 
been used in a variety of cosmographic applications; an overview of past and present work in the 
subject can be found in 0. Lahav's contribution to these proceedings. 

2 Weighted Harmonics and Redshift Distortion 

We start by defining a weighted spherical harmonic decomposition of the flux-limited density field, in 
redshift space as, 



where N g is the number of galaxies in the survey, Y\ m is the usual spherical harmonic, and f(s) is an 
arbitrary radial weighting function. 1 In the special case of constant weighting, one simply obtains the 
projected angular distribution on the sky which is undistorted when the summation in Equation (2) 
is carried out in redshift space; however, in the more general case where f(s) varies with s, peculiar 
velocities will distort the weighted harmonics. 

In order to relate the redshift harmonics given by Equation (2) to their real space counterparts, 
we first rewrite the summation in Equation (2) as a continuous integral over the density fluctuations 
in redshift space, 



In Equation (3), the selection function is assumed to be normalized such that f drr 2 cj)(r) = N g /uj 
where u is the solid angle subtended by the survey. Notice that the selection function in Equation (3) 
is evaluated at the galaxy's distance, not redshift, because if the catalogue is flux limited the probability 
of a galaxy being at redshift, s, will be proportional to the selection function evaluated at the galaxy's 
actual (albeit unknown) distance, i.e. oc (f>(r). Next, we note that, by construction, n s (s)d 3 s = 
n R (r)d 3 r, where n s (s) and n R (r) refer to the densities in redshift and real space and that if the 
perturbations induced by peculiar motions are small, then we can perform a Taylor series expansion of 
all redshift quantities to first order in the density fluctuation, e.g., /(s) ~ f(r) + (U(r) — V 6 s • r), 
where we follow the notation of Kaiser [15] and define U(r) to be the radial component of the peculiar 
velocity field, i.e., U(r) = v(r) • r; V 6 s is the peculiar velocity of the observer at r = (we adopt 
units where H = 1). 

Thus, expansion of Equation (3) yields, 




J2f(si)Yi m (si) , 





(3) 




(4) 



1 In our analysis f(s) is required to be continuous in its first derivative and to vanish at infinity; this simplifies the 
analysis by eliminating surface terms that arise when f(s) has a discontinuous boundary. 
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Abstract 



We examine the nature of galaxy clustering in redshift space using a method based on an 
expansion of the galaxian density field in Spherical Harmonics and linear theory. We derive a 
compact and self-consistent expression for the distortion when applied to flux limited redshift 
surveys. The amplitude of the distortion is controlled by the combination of the density and bias 
parameters, (3 = Q° 6 jb as well as the shape of the real space power spectrum, P(k) (characterized 
by a shape parameter T), and its normalization, <jg; we exploit this fact to derive a maximum 
likelihood estimator for [3, T, and as. We check our formalism using TV-body simulations and 
demonstrate it provides an unbiased estimate of [3 when the amplitude and shape of the galaxy 
power spectrum is known. Application of the technique to the 1.2 Jy IRAS redshift survey yields 
f3 = 0.94±0.17 and T = 0.17±0.05 (l-c) when ag is held fixed at its best value as determined from 
the real space correlation function. Allowing ag to be a free parameter, we find [3 = 0.47 ± 0.25, 
T = 0.15 ± 0.05, and a 8 = 0.81 ± 0.06. 

1 Introduction 

The clustering of galaxies in redshift space appears systematically different from the clustering that 
one would observe in real space. On large scales, the coherent motions of galaxies tend to enhance 
structures tangential to the observer's line-of-sight; conversely, on small scales the peculiar velocities 
smear out structures along the line of sight giving rise to so-called "Fingers-of-God." Consequently, 
galaxy clustering is inherently anisotropic in redshift space. In the context of linear theory, Kaiser [15] 
was able derive a simple expression for the angular dependence of the distortion in terms of the Fourier 
coefficients of the density field, 



where fj, KL is the cosine of the angle between the wave vector k and the observer's line-of-sight. In 
Equation (1) and in the remainder of this paper, we adopt the notation (3 = i7°' 6 /6, where fi and 
b are the current density and bias (assumed here and below to be independent of scale) parameters 





